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O ■ Abstract 

The instantaneous formulations for the relativistic Bethe-Salpeter (BS) and the radiative 



transitions between the bound-states are achieved if the BS kernel is instantaneous. It is shown 
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that the original Salpeter instantaneous equation set up on the BS equation with an instantaneous 

(N ■ 

kernel should be extended to involve the 'small (negative energy) component' of the BS wave 

CO 

functions. As a precise example of the extension for the bound states with one kind of quantum 
number, the way to reduce the novel extended instantaneous equation is presented. How to 
guarantee the gauge invariance for the radiative transitions which is formulated in terms of BS 



wave functions, especially, which is formulated in the instantaneous formulation, is shown. It is 



also shown that to 'guarantee' the gauge invariance for the radiative transitions in instantaneous 
formulation, the novel instantaneous equation for the bound states plays a very important role. 
Prospects on the applications and consequences of the obtained instantaneous formulations are 
outlined. 
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I. INTRUCTION 



It is well-known that Bethe-Salpeter (BS) equation is a relativistic integration equation in 
four dimensional Minkowski space in quantum field theories and it is used to describe bound 
states. It was established by analyzing thepoles and corresponding residues of the relevant 
four- leg Green functions many years ago[l[. Thus BS equation has a solid foundation on 
quantum field theories. The integration kernel of a BS equation (the binding interaction for 
the bound states) corresponding to the relevant Green functions, in principle, contains an 
infinite series of Feynman diagrams, whereas when establishing a precise BS equation, certain 
truncation on the infinite series of the Feynman diagrams must be made by assumption(s) 
i.e. the kernel just corresponds one or a few of Feynman diagrams. Besides computing the 
spectrum of binding systems, furthermore to formulate the transitions between the bound 
states with corresponding BS wave functions was also achieved , so that the applications 
of BS equation have been extended widely. 

Salpeter is the first author, who establishes the precise relation between four- dimensional 
BS equations and three-dimensional Schrodinger equations in describing the bound states 
by means of the so-called instantaneous approximation j^J when the kernel of BS equation 
is 'instantaneous'. In fact, such as positronium, muonium and atoms etc for instance, are 
electromagnetic binding systems, and Coulomb gauge commonly is known to be far superior 
for low momentum exchange to describe the binding of the systems, whereas the leading 
term in the gauge for the kernel is exact instantaneous. Thus if one uses BS equation and 
the kernel is truncated up to leading term only, then they are good examples of the kernel 
being instantaneous exactly. Since there are more mature methods to solve a Schrodinger 
equation than that to solve a four dimensional BS equation, so the so-called instantaneous 
approximation has been adopted widely in various applications of BS equations. Hence in 
literature, when the kernel indeed is instantaneous as realized by Salpeter, quite frequently 
the instantaneous approximation on the BS equation was made and the obtained equation 
with certain 'corrections' to the Schrodinger one generally is called as Salpeter equation. 
Since the original heavy-quark potential model p is on non-relativistic Schrodinger equation 
in three dimension, so people would like to set the potential model on a more solid foundation 
as on a ground of quantum field theory i.e. based on Salpeter achievement to put the model 
on a relevant BS equation with an instantaneous kernel, thus the Salpeter equation has also 
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received additional attentions. 

In fact, heavy-quark potential model itself has made a lot of phenomenological achieve- 
ments since it was proposed. Effective theories, such as heavy quark effective the- 
ory (HQET)jfJ, non-relativistic quantum electromagnetic dynamics (NRQED)[]| and non- 
relativistic quantum chromodynamics (NRQCD)0| etc, were established later than potential 
model, but they have also achieved wide successes too. Whereas in the framework of effective 
theories, when one or more bound states (mesons) are involved in the considered processes 
such as decays and production, they are treated finally to relate to certain local operator(s) 
being sandwiched by the bound states with proper coefficient (s). While the coefficients, be- 
ing of perturbative nature, may be calculated by matching them with those of the underlying 
quantum field theory, e.g., the full QCD in the cases for HQET and NRQCD at a suitable 
high energy scale, whereas the matrix elements i.e. the operator(s) being sandwiched by 
the bound states, being of the non-perturbative nature, can be determined by fitting experi- 
mental data phenomenologically or by lattice QCD etc, such non-perturbative computations 
directly. The non-perterbative computations are difficult and available for limited cases al- 
though they are in progress as the computer power is increasing. Alternatively, sometimes 
phenomenological models may achieve quite accurate results in the estimates of the decays 
and production processes. The heavy-quark potential model, inspired by QCD, is one of 
them and 'quite powerful'. Thus the potential model is still used as a 'tool' to calculate 
decays and production quite often. With it, one not only calculates spectrum of the binding 
system and the 'static' matrix elements but also the matrices for transitions. In this 'di- 
rection', the potential model has been applied and tested quite widely in dealing with the 
non-perturbative nature of the combining constituents (quarks) into bound states (mesons) 
and the transitions. It indeed has obtained quite a lot of satisfied results. Therefore, to put 
the phenomenological potential model on a solid ground of quantum field theory has been 
an interesting topic for quite long time. In fact, in the framework of quantum field theories 
themselves, the problems were solved in many years ago that the bound states are described 
in terms of BS equation and the transitions between the bound states are formulated with 
BS equation solutions (wave functions) so the 'shortest' way is e.g. to embed the 

potential model for heavy quarks into the BS framework on quantum field theories QCD. 

There was an important progress after Refs. . It is about Abel gauge invariance on 
radiative transition formulationP. llfl| . Namely, the so-called irreducible diagrams appearing 
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in the formulation of the transition matrix elements should be determined in certain way 
by the BS equation kernel accordingly: the so-called irreducible diagrams for the transitions, 
which are truncated from a series of the relevant Feynman diagrams, should match to the 
truncation for the kernel of the BS equation exactly. 

For a fundamental process in quantum field theories, the gauge invariance plays a very 
important role, and it is guaranteed perturbatively in Alelian and Non-Abelian cases. In 
non-Abelian cases the gauge invariance is much more complicated when bound states are 
involved in processes. Here as the first step we restrict ourselves to focus the Abelian cases in 
a moment. To guarantee the gauge invariance, although with bound states being involved, 
to work out the precise relation between the irreducible Feynman diagrams for transitions 
and those for BS kernel is not too complicated. In Abelian cases, it is known that for a 
fundamental process the problem is so simple that the gauge invariance is guaranteed, as 
long as all of the Feynman diagrams (either tree or loop ones) in a given order are taken into 
account. But for a process with one or more bound states being involved, to guarantee the 
gauge invariance is not so straightforward, even one may 'trace' them to relevant fundamental 
processes accordingly. It is due to the fact, as indicated by BS equation, that to form a 
bound state means an infinite series of certain Feynman diagrams selected by BS kernel are 
summed, and in the meantime a lot of diagrams out of the selection of BS kernel (there 
are always some in each order) are dropped, therefore, without careful treatment, the gauge 
invariance will be lost. The authors of Ref. II solved the problem for certain cases first, and 

n n 

the authors of Ref. [10] generalized it into general cases. Since Ref. |H| is in Chinese, so we 
will repeat the key points of Ref. briefly in a suitable place in this paper. 

When considering the decays, such as the decays B c — > J/xjj + • • • etc, and there may 
be a great (even relativistic) momentum recoil in the decays due to the mass difference 
m^ c 3> mj/^p (here J/ip is a bound state of a pair of charmed quark c and antiquark c 
and B c is the ground state of c and b guarks), then it is needed to invent a method to 
dictate the recoil effects in he avy -quark potential model properly. To take into account the 
recoil effects properly, in Refs.|lll.ll2|. the so-called generalized instantaneous approximation 
was proposed in terms of BS wave functions and the relevant formulation for the decays. 
The so-called generalized instantaneous approximation was to start with the formulation 
in Refs. for the decays first, which is relativistic so the recoil effects can be treated 
properly no matter how great they are; then it was extended with further approximations to 
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take a generalized instantaneous approximation on the formulation for the transition matrix 
(amplitude) in whole. Finally a formula for the decays was obtained, which may dictate the 
recoil effects well but only Schrdinger wave functions appear . Whereas later on when we 
tried to apply the generalized instantaneous approximation to an electromagnetic transition 
with the achievements, we found the gauge invariance was lost under the approximation. 

In order to solve the fresh problem, i.e., to 'recover' the gauge invariance, we re-examined 
the generalized instantaneous approximation and the original Salpeter equation as well, and 
have found that some un-necessary approximations were made when Salpeter derived the 
equation and when the authors of Refs. jlll. I12I] made the generalized instantaneous approx- 
imation on the transition formulation respectively. We re-derived the equation without the 
further approximations, and a novel instantaneous version of the equation (a group of cou- 
pled equations) to describe the bound states and the formulation for the transitions between 
the bound-states are achieved, as long as the kernel of the BS equation is instantaneous. The 
instantaneous formulation for the transitions is related to the novel instantaneous equation 
directly. In the paper, after some outlining the necessary results briefly, 
we present the re-derivation on the novel instantaneous equations, the 'exact instantaneous 
formulation' for radiative transitions and specially how to guarantee the gauge invariance 
for radiative transitions. Namely we precisely check the gauge invariance of the novel in- 
stantaneous formulation. In the procedure, one may see that although the BS kernel is 
instantaneous as the start point, the so-called 'small (negative energy) components' of the 
instantaneous BS equation (the novel Salpeter equation) are important to guarantee the 
gauge invariance, and cannot be dropped simply. 

This paper is organized as follows: in Section II, we take an example: a fermion and an 
anti-fermion (a mason is the case), briefly to review the BS equation, the four dimensional 
formulation of a radiative transition ^, 3] and to show how to keep gauge invariance for 
themP, Q|. In Section III, we take the same example re-derive Salpeter equation from 
BS equation but with less approximation i.e. to keep all of the components including the 
so-called small components (but not as done by Salpeter to drop the small components). 
Especially, for later usages, we also convert the instantaneous formulation in C.M.S. of 
the bound state into a 'covariant' instantaneous formulation. Then we present the brief 
derivation and final results on instantaneous formulation for the radiative transitions between 
bound states from the four dimensional one to effective three dimensional one. We also check 
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the gauge invariance for the instantaneous formulation in terms of the four dimensional one. 
In Section IV, we discuss the obtained formulations and their consequences for effective 
theories etc when the BS equation has an instantaneous kernel. We make prospects on the 
applications of the instantaneous formulations. In Appendices we present some formulae 
used in the derivation. 

II. FOUR DIMENSIONAL BS-FORMULATION FOR TRANSITIONS BETWEEN 
BOUND STATES AND ITS GAUGE INVARIANCE 

In this section, we show the four dimensional formulation for covariant BS equation and 
possible radiative transitions between the bound states in terms of the BS functions and 
show how to guarantee the gauge invariance for it etc. 

A. The BS Equation 

In general, for a binding system of a fermion and an anti-fermion, the BS wave functions 
Xp^(xi,x 2 ) and Xp^( x 1i x 2) are defined 

Xp&i,x 2 ) =< 0|T[V(zi)^Or 2 )]|P£ >= e wx XPi; {x) , 
Xp&i,x 2 ) =< P$\T[$(x 2 )il>{x 1 )\\Q >= e- wx xp&) , (1) 
X = aiXl + a 2 x 2 , x = xx-x 2 , «i = , « 2 = , 

where X, P are the CM. coordinator and the four momentum of the bound state; 
m 1 ,m 2 ,Xi,x 2 are the masses and coordinates of the 1-quark (q{) and the 2-antiquark (q 2 ) 
respectively. The wave functions in momentum space 

^;Xp((q) = I d*xe- iqx xpz(x), ^=xps(?) = / d 4 xe iqx xpd x ) , 

P = Pi + P2 , q = OL2V1 - "1^2 , (2) 

where q, u are the relative four momentum in the bound state and the energy of the bound 
state respectively. In Eqs. ([l|2|) £ denotes the possible and 'inner' degrees of the freedom 
(quantum numbers) of the bound state (such as the mesons A and B etc). Since £ does not 
affect our derivations and theorems in the paper, so we drop it out in all of the statements 
later on. 
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The BS equation generally may be written in momentum space as [If 
J d 4 q0(P, q', q) - V (P, q' , q)} X p(q) = , 

/ A'xp(g') (0(P, q', q) ~ V(P, q', q)} = . (3) 

The 'differential' operators O and O here acting on the BS wave function under integration 
accordingly are defined as the follows: 

J d i q'd(P,p,q)xp(q) = J <Pq5\p - q)[Sf {a x P + p)]- x Xp(ff)[5f (-a 3 P + q)}' 1 , 

J d*pXp(j>)*0(P,p,q) = J d^p-q^Sfi-^P + q^XPmsfi^P + p)]- 1 . (4) 

To be explicit, here we use Sf with the superscripts "(1)" and "(2)" to denote the propaga- 
tors of the quark qi and the anti-quark g 2 respectively. Throughout the paper even without 
explicitly explaining, we will use the superscripts to denote the relevant quantities or opera- 
tors, i.e., to denote those for the quark qi and the anti-quark q 2 respectively as here. While 
V(P,q',q) and V(P,q'q) are binding interaction i.e. the BS kernel. In addition, the BS 
wave function Xp{q) satisfies the normalization condition: 



d - 3 ^~Tr \x P (q)^ ([^7 V)] (1) [Sj 2) (-P2)]- V(g - q') + V(p, q, q')) X M) \ = 2iP , 

(5) 



(2vr) 4 r^aPo 



here Pq = uj is the zero component of P. 

If the 5 function in Eq.(J3J), contained in operator O (P,p, q) as indicated in Eq,(JD), is just 
integrated out, the Bethe-Salpeter (BS) equation corresponding a quark- ant iquark bound 
state can be written as: 

r d 4 k 

(fa - mi)x P (g)(/ 2 + m 2 ) = J j^Xa V { p , 9> k )Xp(k) , (6) 

or 

Xp(q) = -r^— I -jP^ViP, q, k) XP (k)—^- , (7) 
fx — m\ J (27r) 4 p 2 + m 2 

which is independent on the specific reference frame because it is Lorentz invariant. 



B. The Formulation for Radiative Transitions 

As an example of the radiative transitions, let us consider a decay A(P) — > B(P') + j(Q), 
and we will take this example to show the formulation for the decays in terms of the BS 
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Pi 7p I p[ 




m% = r«2 



FIG. 1: Feynman diagram corresponding to the quark in the 'meson' A(B) emitting the photon 
for the transition A(P) -> B{P') + j(Q). 

Pi = p'i 




FIG. 2: Feynman diagram corresponding to the anti-quark in the 'meson' A(B) emitting the photon 
for the transition A(P) -> B{P') + j(Q). 

wave functions and how the electromagnetic (an Abelien) gauge invariance is guaranteed 
throughout the paper. Here the 'meson' A (a bound state) consists of a quark q\ : p\,ra% 
and an anti-quark q 2 : p 2 , rri2, and the 'meson' B (another bound state), similar to A, consists 
of a quark q[ : p'\,m'i and an anti-quark q' 2 : p' 2 ,m' 2 . To be responsible for the decay, the 
Feynman diagrams for the lowest order are described as Fig^and FigEl 

Since the concerned process is an electromagnetic transition, and if it is 'automatic' 
decay, the meson A must be an excitation state of the meson B, i.e., we must have mi = 
m' 1 ,m2 = m' 2 . For more specific, we further assume that in the mesons, the quarks qi,q[ : 
(pi, mi), (pi, m'i), carry charge Qi in the unit e and the anti-quarks q 2 , q' 2 '■ (P2, (j>2, m'2), 
carry charge Q 2 in the unit e. In fact, in many discussions, proving the gauge invariance 
and the other derivations, the contributions from quark qi Fig^and those from anti-quark 
qi FigOare independent, thus to investigate the decay amplitude, either the quark emission 
of a photon FigQor the anti-quark emission FigOmay be taken as 'a representation' of the 
transition (to consider one of them is enough, the other is almost the same just to repeat 
once more if one consider both of them for completeness, so here we say 'a representation'). 
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In Refs. DD , how to write down the S-matrix element for this kind of decays in terms of 
BS wave function is given: 

< p', q e\s\p >= jf^ 6 8\p' + q- py. < py^p >, 

here is the polarization vector of the photon and the matrix element for the current 
(Figs.1,2): 

< P'\J,\P >= / d 4 q'J d 4 qTr{xpiq') [Gf{P'q'; Q; P, q) + G®(P'q'; Q; P, q)] X p 

= jhfl A 7 ^ {xA^Q^XpiQ^Sfiq' ~ a' 2 P')}- 1 ■ ~ 4 + ol % P> - a 2 P) 
+xAq')[Sf\a l P + q)]- l Xp{q)Q2l^ • 5\q - q' + a,P - a[P>)} , (8) 

namely corresponding to the irreducible Feynman diagrams Figs. (jll2j) . the operators (fac- 
tors) Gf{P'q'] Q; P, q), (i = 1, 2) are defined byQ 

xAq , )G^(P'q';Q;P,q) X p(q) = xAq')Qu^Xp(q)[Sf\q' -c^P')]"^^?' '+4? -a 2 P) ; 

xAq')G^(P'q';Q;P,q)xp(q) = xAq')[Sf ] {^P +q)]- 1 x P {q)Q 2l ^ -S^q-q' +a 1 P-a' 1 P') . 

Here the quarks' momenta, p\,p\, and anti-quarks' momenta, P2,P 2 , are related to the total 
momenta P, P' and the relative ones q, q' of the initial meson A and the finial meson B 
respectively: 

r> , _ m x _ m 2 

p 1 = a 1 P + q, a>i = , p 2 = a 2 P - q, a 2 = 



mi + m 2 f^i + wi2 

p = a p + q a 1 = —j- T , p 2 = a 2 P - q , a 2 = —j- T . (9) 

m\ + m 2 m\ + m 2 

In Eq.(JHj) 'Tr' denotes the trace for the 7-matrices is taken. Xp{q) an d Xp'(<z') are t ne 
relevant BS wave functions of the mesons A and B; Sf(i = 1,2) are the propagator of the 
quark and anti-quark. Q = P — P' is the photon's four momentum. Since the decay is 
electromagnetic, i.e., the flavor (s) is(are) conserved, so we have m\ = m'i and m 2 = m 2 , 
that «i = a[, a 2 = a 2 and a[P' — a>iP = —aiQ, {i = 1, 2). Thus later on we simplify Eq.(jHJ) 
by taking rru = and m 2 = m 2 always. 

In Refs.Q, Q, the operators G^{P'q';Q; P,q) and Gf>{P'q';Q;P,q) are obtained by 
analyzing the Feynman diagrams and to be sure only the 'irreducible' diagrams are involved. 
In fact, the amplitude being gauge invariant also means the operators have very definite 
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relation to the BS equation (the 'differential operator 0(P,q,q') or 0(P,q,q'Y and the 
kernel V(P, q, q')). Let us show the relation of the operators to the BS equation in the 
following subsection. 

C. The Electromagnetic Gauge Invariance of the Radiative Transitions 



According to Refs.|9lllQl|. when a transition current matrix element, in which there is(are) 
bound state(s) being involved, is formulated in terms of BS equation, such as Eq.(jHJ), and 
if one requests it is gauge invariant (the current is conserved), then the transition opera- 
tors G^(P',q';Q;P^a)(i = 1,2), which appear in Eq.fjHJ) and correspond to the irreducible 
Feynman diagrams[l7], must 'match' with the BS equation accordingly in an exact way [18]. 
Therefore, to determine a specific G^(P',q';Q; P,q) and to guarantee the matrix element 
being gauge invariant, we need to explore the match rule precisely. 

Now let us 'explore' the match rule and show the specific Gff(P', q'\ Q; P, q) happens to 
be satisfied with the match rule, when the BS wave functions adopted in Eq.fjHJ) correspond 
to the solutions of the equation with the 'lowest' order kernel. 

To explore the 'match rule', first of all, let us introduce an external classical field A e * and 
in the underlying theory add a coupling term: 

J^(x)A^(x), (10) 

to the Lagrangian, here Jj? m is the electromagnetic current of the field theory. Then being 
parallel the case without the A e *, one now may establish a BS equation with the external 
field Aff through the coupling term Eq. (PTj|) . and the BS equation operators A< ^ and V A ^ 
may be obtained accordingly. With the preparations, the 'match rule' can be stated as that 
G M , the irreducible Feynman diagrams, should be exactly obtained as the the following way: 



{G»}= E Gf{P',<j-Q-P,q) = 6{ ° ] J - } U^o. (11) 

i=l,2 0A fi 

The Eq. (JllJ) . in fact, can be explained as that the irreducible Feynman diagrams, G^, should 
be obtained by the way that the photon line attaches onto each charged line in the BS 
equation operators 0(P,q',q) — V(P,q',q) once in turn. To be precise, let us take the 
example to explain Eq. ljllj) and its meaning. 
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With the match rule between G fl (P',q';Q;P,q) and BS equation, let us turn back to 
examine the specific Gf{P', q'\ Q; P, q) in Eq.flHJ). If the BS kernel V(P, q', q) for the bound 
states P, P' of quark and anti-quark is just one gluon exchange 

d 4 q'xp(q')V(P,q',q) = J A'7? XpW)^^ ~ q) 



or 



J d 4 qV(P, q', q) XP {q) = J d 4 q^ XP (q)^ D f (q' - q) , (12) 

i.e. Df is the gluon propagator, and in the underlying theory the quark carries charge Qi 
in unit of e and the anti-quark carries Q 2 that means in the electromagnetic current J£ m 
Eq. ([T0|) there are quark current QiJ^ and anti-quark Q 2 J^\ hence by Eq. ffTTj) for the 
transition the irreducible Feynman diagrams G^(P',q';Q; P,q) is generated by attaching 
the photon line onto the quark line in 0(P, q', q) — V(P, q', q) (Fig. la) and G^'(P', q'\ Q; P, q) 
is generated by attaching the photon line onto the anti-quark line in 0(P, q', q) — V(P, q', q) 
(Fig. lb) as indicated in Eq.fjHJ). Note here that since V(P, q', q) is neutral so according to the 
match rule there is no contribution from V(P, q', q) at all, thus the photon line can attach 
onto those lines in Q(Pjj', q) of the BS equation only. 

Now following Refs.pl, llflj . let us show the gauge invariance, i.e. the current matrix 
element is conserved. 

To apply the equations Q = P — P' and the general identity 



- iCtl® = [Sfioc.P + q)]- 1 - [S { ;\a[P' + qT 1 , 
Wlf = [Sf(-a 2 P + q)]- 1 - \Sf\-a' 2 P' + q')}' 1 , (13) 



we have 



-iQ" < P'\J,\P >= I d 4 q' I d A qTr {xP'{q')Q\ {[Sf^P + q)]' 1 
-[Sf\a[P> + q')]- 1 ) X p(q)[Sf\-o/ 2 P> + q')}- l 5\q - q> - a 2 Q)- 
xAq')[Sf\ ai P + q)Y l Xp(q)Q2 ([Sf(-a 2 P + q)]' 1 - [Sf{-a> 2 P> + q')]^) 5\q - q' + Ul C 
= J dS I d A qTrxp>{q') {Qi [ct(P, q' + a 2 Q, q) - fr(P', q>, q - a 2 Q)] + Q 2 [■ ■ ■]} X p(q) ■ 

Considering the property of the BS kernel V indicated by Eq. (jl2j) 

XP>(q')\?(P, q' + a 2 Q, q) XP (q) = xAq')^{P\ </, q ~ a 2 Q) X p(q) , 
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and 



xp'W)v(p, q - aiQ, g)xp(q) = XP'{q')^{P', q, q + &iQ)xp{q) , 



(15) 



then Eq. (jl4j) can be rewritten 



-iQ^ < P'\J^\P > = 
J dS I d 4 qTr {QiXp'(?') [#(P, q' + a 2 Q, q) - t(P, q> + a 2 Q, q)] X p(q) 



- QiXP'(q') p(P', q', q - a 2 Q) - V(P', q', q ~ a 2 Q)} X p(q) 
+Q2Xpiq') q' - UiQ, q) - V(P, q' - ai Q, q)] XP (q) 
- QiXfW) p{P', </, q + <*xQ) - ^(P', q', q + x P (q)} 



(16) 



It is easy to see when the BS equations Eq.® are applied to Eq. (fT(3|) for Xp'(p) an d Xp{q) 
respectively, the matrix element of the current is conserved 



i.e. the Abel gauge invariance is proved. 

Since the charge values Qi and Q 2 of the quark and anti-quark are independent, and one 
may also see from the above that the terms proportional to Q\ and Q 2 are independent too, 
thus to prove the gauge invariance, it is enough to show one case, for instance, for the terms 
which are proportional to Q\ only. Namely it is enough to show the terms proportional to 
Qi themselves are gauge invariance (for the terms which are proportional to Q 2 to show 
their gauge invariance is similar). Later on to shorten the presentation of the derivation, 
we only present those terms either proportional to Qi or to Q 2 precisely, and just add the 
statement as 'the other terms are similar' for instance. 

III. THE INSTANTANEOUS FORMULATION FOR RADIATIVE TRANSI- 
TIONS BETWEEN BOUND STATES AND ITS GAUGE INVARIANCE 

Before deriving the instantaneous formulation for the transitions between bound states, 
first of all, we need to do some preparations, namely we need to re-write the 'instantaneous' 
BS equation in covariant formulation in a reference frame, in which the bound state may be 
moving. 



-iQ» < P'\J tl \P>=0, 



(17) 
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A. The BS Equation with an Instantaneous Kernel and Its reduced Equation 



The so-called 'instantaneous approximation' was realized by Salpeter first, and it is not 



Lorentz covariant 



4] 



In his original paper and in most literature, indeed an approximation 
was made in the rest frame of the concerned bound state. It says when the BS has an 
instantaneous kernel i.e. the kernel of BS does not depend on the time component of the 
relative momentum k , such as to take a simple form 

V(P", k^q^ =o ~V{\k - ? |), (18) 

then as a main result of the 'instantaneous approximation', the four dimension BS equation 
can be reduced into a three-dimension Schrodinger equation jj]. 

Whereas here when BS equation really has an instantaneous kernel, we reduce the BS 
equation without approximation in certain sense. Especially, we are considering the radiative 
transitions between the bound states, and we find that with an approximation as done by 
Salpeter the electromagnetic gauge invariance for the transitions cannot be guaranteed, 
thus we need to start the reducing of the BS equation with an instantaneous kernel from 
beginning, but to reduce the 'instantaneous equation and the wave function' in such a way 
that is the so-called instantaneous but Lorentz covariant. 

To do so, first of all we divide the relative momentum q into two components, q.. and q ± , 
a parallel component and an orthogonal one to the bound state momentum P (P 2 = M 2 ) 
respectively: 



9" = 



where = ^rP M and qp ± = g M — qp . Correspondingly, we have two Lorentz invariant 
variables: 

_ (P ■ q) _ r- 2 - r - r 

Qp- M » QPr-yQp-Q -y-Qp ± - 

In the rest frame of the concerned bound state, i.e., P— 0, they turn to the usual components 
qo and |"(f|, respectively. Now the integration element of the relative momentum k can be 
written in an invariant form: 

d 4 k = dkpkp T dkp T dsd(p, (19) 

where is the azimuthal angle, s = ■ The interaction kernel Eq. (jl8j) can be 

re-written as: 

V(P^k^q») = V(\k^-q» P± \). (20) 
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Similar to Ref. 41. we introduce the 'instantaneous' BS wave function: 

M^^'/^.i), (21) 

and integration 

= / %ly dS v (K ± - & ± \)M*rpJ , (22) 

then the BS equation EqEI becomes 

Xp(Q P{l ,qpJ = SfWMqpJSfi-pt), (23) 

where S^\pi) and (— p 2 ) are the propagators of the quark and anti-quark respectively. 
They can be decomposed as: 



- iJSfiJpt) = + , (24) 

3 Jqp + OjM — u iP + ze Jqp + otiM + uj iP — le 



with 

UiP = \/™l + (lp T , A± (q P± ) = — 



p 



M 



(25) 



where J = 1 for the quark(i = 1) and J = —1 for the anti-quark(i = 2). Here Af p (qp ± ) 
satisfies the relations: 

A£(#J + A£(g£J = ^ , A±(^ L )^A±(^ L )=A±(^J, A±(^J^AS,(rfJ=0. 

(26) 

Namely may be considered as projection operators (the energy (P)-projection operators), 
while in the rest frame of the bound state, they turn to the energy projection operator. 
For later discussions let us define fp^iq^) as: 

#(#J = Atp(Q P J^p(<lpJ^fp(<lpJ , (27) 
then for the BS wave functions <fip{qp ) we have: 



Vp(9p, ) = ( PV(Qp l ) + Vp (Qj>, ) + ( Pp + (Qp i ) + <^p («p. 



Since the BS equation kernel is instantaneous, we may integrate out qp (contour integra- 
tion) on both sides of Eq.(|25|). and obtain 

f „ , *tp(<£Jvp(<£Jt4p(<&J ^ipWpJvpWpJ^pWpJ foa , 

( fP{ ( lp ± - 7TT v 771— ; v , l^oj 

± (M — ujip - uj 2 p) {M + uip + u)2P) 
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and applying the energy (P)-projection operators Af P (q P ) to Eq.(J3HJ) further, we obtain 
the coupled equations: 

( M - uj 1p - u 2P )^ P + {q P J = AUlP ± )Vp(q"pJ A tp(qk) > ( 29 ) 

(M + u 1P + u 2P )<p p -(q P J = -ATMJvp(qpJ^p(QpJ > ( 3 °) 

<P + p-(tiJ = Vp + (tiJ = 0- (31) 

The normalization condition Eq.(jSJ) now terns to the follows accordingly by the positive 
and negative functions as: 



2P . (32) 



qldq T _ ++ j) ++ i> j) 

-f —tr 09 — 09^ - 09 — 09 — 

(2tt) 2 [_ M M M M 

Note that Salpeter and many authors in literature considered that, of the above coupled 
equations Eqs. ()29l30l31|) . the equation Eq.([29j) is 'dominant', owing to having a very small 
coefficient (M — ui P — dj 2 p) (weak binding) on the left hand side, and it is 'projected' by 
positive-positive energy project operators (Af P (q P , ) and Ajp(gp ± )) (it is why people call it 
as positive energy equation). Thus the authors highlighted on it only and dropped the rest 
equations out at all. With the approximation, they showed that if expending the equation 



Eq.flU to 

accordingly 



;he order of -rW further, the equation may be turned into a Schrodinger one 



r lo|. Whereas, as emphasized in Introduction, we realized that if only keeping 
Eq.()29J) and its solutions as wave functions to be applied to the radiative transitions, the 
gauge invariance will be violated. Thus here in the next two subsections we will show the 
formulation of the transitions between bound states in four dimension how to turns into an 
instantaneous one and accordingly how to keep the gauge invariance when the formulation in 
four dimension is gauge invariant precisely. It is just by considering all the coupled equations 
Eqs. f|29l3QI3H) 'equally', and one will see that all of the equations play important roles in 
guaranteeing the gauge invariance for the radiative transitions between bound states. 

To explore the contents of the Eqs. (j29l3UI31|) . here let us take a heavy quarkonium (or 
positronium) as samples to establish the coupled equations for the components in the in- 
stantaneous BS wave functions <p P {q P± ) and to show the equations are consistent. 

Since a heavy quarkonium (QQ) is of equal masses (mi = m 2 , i.e., a± = a 2 = |), so its 
physical states have definite quantum number of J PC , thus its equations should be given 
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according to the quantum number precisely. Here we take the state with J PC = 0~ + (total 
spin 5 = 0, and total orbital angle moment L = 0) as an example to show the main feature 
of the equations Eqs. (|29l3QI31|) . Of the components, the general relativistic wave function 
for the bound state with the quantum number J PC = h (in the center of mass system) 
can be written as: 



7o¥>i(9) + f2(l) + 4^0) + 7o(/±<^4(<?)] 7s , 



(33) 



where q± = (0, <?). The equation Eq. (f3"T]) . in fact it acts as constraints, give the strong 
conditions precisely: 

<pi(q) = -mnpi(Q), <Ps(q) = • (34) 



Then we can rewrite the relativistic wave function of state + as: 



iofi(i) + <f2{ ( i) -7o^± — <pi(q] 

mi 



7s- 



(35) 



To apply the definition Eq. (j2"Tj) to Eq. (j33|) . the corresponding positive and negative projected 
wave functions may be written as: 

1 



2<jj\m\ 



WiViiQ) +m 1 (p 2 (q)) (7 mi + coi -7o(/j_)75 



and 



fi So (Q) = — — — ((ujipi(q) -mi¥> 2 (9)) (7o™i - o>i -Jad±)j5- 



(36) 



(37) 



From the reduced BS, Eqs.(J2Hl EDJ with a QCD inspired instantaneous kernel (relating to 
V s and V v ) for the rj in the equations and by taking different traces for the F matrix, finally 
we obtain the couple equations for the components <fii(q) an d ^>z{9)' 



-/ 



dk 
(2tt)3 



(M — 2u)i)u 1 Wi<pi(?) +mi<p 2 (q) 
m\(V s + 2V V ) - q ■ kV s ) yi(h) + uj 1 m 1 (V s - 4V v )(p 2 {k) 



(38) 



-/ 



dk 



m((V s + 2V V ) - q ■ kV s ) ^i(fc) - u imi (V s - AV v )tp 2 (k) 



(39) 



i.e. coupled equations for the eigenvalue M and eigenfunctions, the wave functions <pi(Q) 
and ip 2 (q). Note that we will derive and solve the equations with kernels V s (corresponding 
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to scaler confinement) and V v (corresponding to one-gluon exchange) precisely elsewhere 1 131]. 
but here we only show the fact that when all of the equations Eqs. (|29l l3*Ul l3*T|) are taken into 
account, consist coupled equations can be obtained finally without approximation. 

By contrary, for Salpeter equation only the positive energy component equation Eq. (J29j) is 
kept, but further assumption on the spin structure (in four components) of the wave function 
ipiSoiQ) i s niade|4|, that is different from here, as shown, the spin structure is determined 
by the equations Eqs. (|2*9l EM l3*Tj) simultaneously. 

In fact, the above coupled equations are written in the reference frame of the bound 
state itself. Whereas, to give the 'instantaneous' formulation for the transitions between 
the bound states, we also need to have the precise formulation for the BS equation (with 
'instantaneous' kernel) in a reference frame of another bound state. It is similar to do as 
the above, so the precise formulation is put into Appendix A. 



B. The Reduced Formulation for Radiative Transitions between the Bound States 
with Instantaneous Kernel 

The transition matrix element for a process of a bound state emitting a photon can be 
written by the BS wave function as: 

< p'\ j„\p >=< py^p >! + < p'\ j„\p > 2 , (40) 

where 

< P'\J,\P >! = J ^0^Tr {XP>Q11»XPS 2 \-P2)5\P2-P 2 )} , 

<P'\J,\P> 2 = J d -0^Tr{ XP ,S^\ Pl ) XP Q 2lll 5\ Pl - P ' l )} , (41) 

here we have used subscript '1' to denote the photon emitted only by the quark and '2' to 
denote the photon emitted by the anti-quark. In fact, when we consider the gauge invariance 
of the process, these two parts are independent and each one can be treated independently. 
So below when we show how to keep gauge invariance, we pay more efforts on the first one, 
and just to say similarly on the second one for shortening. 

We have made instantaneous approach to the BS wave function and BS equation, fur- 
thermore in this subsection we do similar 'job' for the transition matrix elements. 
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As for the matrix element 
f d^qd^q' 

< P '\J^\P >1= Ql J (2tt)4 Tl " { S '^-P2)v'p'S[{p l )l^S l (p l )T]pS 2 (-p 2 )S 2 1 {-p2)^{P2 -P 2 

= Qi j ' ~^^~ Tr {v , p'S' 1 (p' 1 h^S 1 (p 1 )ripS 2 (-p2)5 A (p2 - p' 2 )} , 

here the integral element d A q = d 3 q p ±dqp and the ^-function can be re-written by the 
relation: 

1 



$\P2 -P 2 ) = -S 3 (P2 P ±-P 2P± ) x — ( 



a 2 P P — q P — a' 2 P' P + q' p - ie 
then with Eqs.(|24 | 125 )1 . the matrix element turns into 

d A qd A q' 



2m ^a 2 Pp — qp — a' 2 P'p + q' P + ie 
' N (12) 



< p '\M p >i= -Qi J fr y Tr {n'p>s'i(p'ih»Si(p 1 )r)pS 2 (-P2)s ;i (p2P± -p' 2PA _) 



(2tt)< 

1,1 1 



x 



2m a 2 P P — qp — a^P'p + q' P + ie a 2 Pp — qp — a' 2 P' ' p + q' P — ie 



.dVV J-, f a+ Ar 

-<3i / /o_^ Tr S v'p'SAp^ I , ; nj — — + 



(2vr) 4 [ iV u "* Xqp + aiM-uJx + ie q P + a x M + u x - ie / 

{ A t A 2 \ 3 , 

x Vp ; T7 -7- + : ti-. o {p 2P± - p 2P± ) 

\ -q P + a 2 M -tu 2 + ie —q P + a 2 M + tu 2 - ie J 

i(i i M 

X- 



2m \a 2 M — qp — a' 2 E' + q' P + ie a 2 M — qp — a' 2 E' + q' P — ie J J 
To integrate over qp and in terms of Eqs. (j29[ IHU|) for the initial bound state, the above 
matrix element becomes the follows 

d 3 q P ±d 4 q' 3 , , ( - / ¥?p + (9p±) 



< P'\J,\P > 1= -Q, j ^ q S 3 (p 2P± -p' 2p± )Tr v'p'SMH 



+ 



(2vr) 4 ^ 2F±J y""-L\ruwy M — a 2 E' — lui + q' P + ie 

<fip + {qp±) a iVp^2 



—a 2 E' + u 2 + q'p - ie (M - lu 1 + q' P - a 2 E' + ie)(a 2 E' + u 2 - q' P - ie) 
AiTipAt (p P ~(qp±) 



(M + uji + q' P - - ie)(a 2 E' — u 2 — q' P + ie) —a 2 E' - u 2 + q' P + ie 

M - a 2 £' + wi + g P - z'ey J ' 
To apply the formulae in appendix A and to carry out the the cantor integration over q' P , 

we have 

Vp/A^p/ {q'pi_)l^V{<lP±) 



< P'\J,\P > 1= Qi y gg^k ^^-^jrr | ^ 



E' — oj[ — lo 2 
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ri'p>Kp>Wp±h» A iVP^~2 



M — u)\ — u)' x — E< (M — lu 1 — lu' — E')(E' +co[ + lu 2 ) 

v' p'K + p'(qp±hn A iVpM v'p'^ P 'Wp±)i^p~((1p±) 



{M + uj\ + uj' - E'){E' - lu[ - lu 2 ) M + ivi + u[ - E< 

Vp,A / 1 p,(g^ ± )7^p~(g P± ) \| 
E' + uj[ + uj 2 j J ' 

Here we have 5 3 (p 2P± - p' 2P± ), u' 2 = u> 2 , A-' 2 ± p/ (q' P ±) = Af P (q P± ), the identity 1 = £J = 
[A! 2 + A' 2 ~) and the equations Eqs. (p29l IHUj) for the final bound state, to apply all of them 
to the matrix element, we finally obtain: 



<P'\J^\P> 1 =Q l J 



d- i q P± d- i g' P± 
(2vr) 3 



5 3 (P2 P ± -p' 2P± )Tr 



L 

M 



1ft p, l^p + 



M + u)i + u[ - E> 



Kp'Wp±)v'p'^p'Wp±h^p (qp±) A ' 2 + p'(<i'p±)v'p'Kp'Wp±h^p + (Qp±) 



M + u x + u{ - E' 

Lp'p7'j fl Afr] P A 2 



M - ux - J x - E' 



If we define: 



M - u x - u[ - E' 



(f' P , i^p 



if, 



ip 



Ai (qp±)ypAt(qp±) 

M + LU X + lu[ - E' 



(43) 



(44) 



K~Aq'p±)v'p'K + Aq'i 



p±j 



M + + — E' 



(45) 



ip' lP , 



a ' 2 + p'(q'p±)v'p'Kp' 



M - L0 t - lo[ 



E' 



(46) 



i>tp = 



Aj(q P ±)r] P A 2 (q PA _) 
M -u-i-u'x-E' 



the amplitude becomes 
<P'\J,\P> 1 =Q 1 



d 3 q P ±d 3 q P± 3 , 

{-qp± - ol 2 P P ± + q P± )Tr 



(47) 



+ip' P ?(q' P± ht,ip P + (qp±) - ip' P ?(q' P ±h^ P ~(qp±) - i> ' P ~ (q p±h^p + (QP±) 

+ftpT {q' PA _)l^p~ (Qp±) ~ <f'p~ (O 1 'px)1i* ( Pp~ (Qp±)] } • ( 48 ) 

To integrate over three dimensional q P± further, then we obtain the instantaneous transition 
matrix element < P'\J^\P >i in three dimensional formulation. To do the same derivation 
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for < P'\J U \P >2, the whole instantaneous transition matrix element becomes 



< p\j^\p >=< py^p >! + < py^p > 2 

= / ~{2^ Tr {^M [ <fif 'r( qpl + a 2^'p±)7/^p + (9p±) + 1 P'p^( ( 1p± + «2-P ' pih^Tp (qp±) 

- Tp'ip>{qp± + at 2 P'p±h»<Pp~(qp±) ~ 4>'tp>( ( lP± + a2P'p±hn¥p + ( ( lP±) 

+ <-p' P ^(qp± + « 2 -P / p±)7mV'ip (qpx) - v'p~(qp± + a 2 P'p±) r r^ P '(qp±) 
+ Q2 v'p^(qp± - aiP'p±)jj¥p + (qp±) + ftV(Qp± ~ a i p 'p±)M4>2p(qp±) 

- *P'tp>(qp± - a i p 'p±)ji i Pp~(qp±) - iP'2P>(qp± - aiP'p±)§^^t + (qp±) 

+ <p'p~(qp± - &iP'p±)j<*p2p(qp±) - v'p^iqp^ - "i^'p±)§<^ P ~(gp±)] 7^} • (49) 

Here similar to the Eqs. (|44M4"7j) . we have defined ibfp and ip 2 pf as follows: 

_ ^t(qp±)vp A 2(qp±) 



lb 



M + uj 2 + uj' 2 - E' 



(50) 



^2P' 



K + Aqp±Wp'K~Aqi 



p±j 



M + lu 2 + lu' 2 - E' 



(51) 



/-+ A > 2P , (qpx)v'p'A'iP'(qp±, 



ip 2 p> 



M - uj 2 - lu' 2 - E' 



(52) 



Ai (qp±)vp A t (qp±) 



2P 



M — LU 2 — lu' 2 — E' 

If one takes the approximation <fp~(qp±) — 0, the matrix element becomes 



< py^p >=< py„\p >! + < py^p > 2 

—^ Tr W 1 ]^ WP'iqPl + a 2 P' P ±)-f^ + (qp±) 



(2n 



+ (p'pt(q P ± + a 2 P' p±h^ip (Qp±) ~ ^'ip'(qp± + a 2 P'p±h^ ( Pp + (qp±, 



- P - ++ P 

<p'p> (qp± ~ aiP'p±)jj<Pp + (qp±) + ( P'p^(qp± - oiiP'px)—^p~(qp±) 



+ P 
~ ^'2~i>{qp± - a 1 P'p±)—(pj, + (q P ±] 



7m 



(53) 



(54) 



Note that even when the approximation (p P (qp±) — is made, from Eg. (j54j) one may 
see that it is still not correct to calculate the matrix element by means of 



< pyjp >=< pyjp >! + < pyjp > 2 
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(2tt) 3 

+ q 2 



r 



<P' P > (qp± - ®iP'p±)jj¥p + (<lP±) 



(55) 



C. The Gauge Invariance of Instantaneous Transition Matrix Element 



In Sec. II, we showed that the transition matrix element is gauge invariant no matter the 
instantaneous approach is applied to, i.e., 

1 



- iQ» < P'\JJP >i- 



(27T) 



A' J d'qTr {xAq')Qx {[Sfia.P + q)}- 1 
\Sf\a\P* + q')]- 1 ) XP (q)[Sf(-a' 2 P' + qT'^q ~ <t ~ a 2 Q)} = . 



(56) 



In this subsection, we show that the instantaneous transition matrix element Eq. (|49|) is also 
gauge invariant, in terms of the gauge invariance in four dimensional formulation Eq. (|56l 
and the instantaneous equations Eqs. (|29l3UI31|) . 

First of all, let us examine how the equation Eq. (J56j) turns to a three dimensional one 
when the kernel of the considered BS equation is instantaneous. To substitute the BS wave 
functions Xp(q) an d Xp'W) with BS equation Eq.(j2HJ), and similar to the above subsections, 
to deal with 5 functions and then to integrate out qp and q' P in the equation by contour 
integral, one immediately obtain 

-iQ» < P'\J,\P >i= / d ~0f T r {s7\-v'z)VP.sY\p' l )r ] p - fj P 'Sf\ Pl ) V pSf (-p 2 )} 

x^ 3 (P2p± -p' 2P± )5(a 2 M -q P - a 2 P' P + q' P ) 

. f d 3 q ± d 3 q' ± 5 3 (p 2± -p' 2± ) r , ++ _ ++ , - —\ 

= -2tt? J j—-^ ^Tr [tpfrriP ~ Vp> Vp ~ Vp'VpJ + Vp"Pp \ • (57) 

Note here that with Eqs. (j29t l3"U |) . we have the part in Tr of Eq. (j57j) : 

{fi'p> + VP ~ Vp^Vp - VP ,( fpt + Vp'¥p~} 

\ A' 2 + p ,(q' P± )rj' p,A'£>,(q P± ) Kp l (^p ± )rj' pi A'{ pl (g' P± ) 

~Vp H ^r, ; — ; Vp 



E' — uj 2 — u[ 



E' + u 2 + u' x 



-, A-ip{qp±)VpA-2p(<lP±) -, A 1P (q PX )r]pAt P ,(qpx) 
'V P' U ~ Tj ^ P' ' 



(58) 



M — uj 2 — uj[ r M + to 2 + oj'i J 

From Eq. (j58|) . one cannot be sure the difference between the second and the fourth terms 
will be smaller than the difference between the first and the third terms, although each of the 
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first and the third terms itself is big and each of the second and the fourth terms is small. 
Practically, we have taken several examples for heavy quarkonium, and indeed we have 
found when <p'p,~ and (p]>~ in Eq.(|57|) are neglected, then we cannot obtain satisfied results 
with gauge invariance. Therefore, when one would like to have the gauge invariance for the 
transition matrix elements, one cannot ignore any of the 'negative energy' components in 
Eq.dH- 

By means of the relations of projection operators Eq. ()26|) . the above equation Eq. (l56l) 
becomes: 



-iQ- < PV„\P > 1= -2*i / dVA1 ffi~ ri±) 7> {^{jW + AD, P (AJ + A 



Z 

M 







Furthermore, with Eqs. 
-iQ" < P'\J tl \P> 1 



, we finally have: 



(27T) 



J d 3 q ± Tr 



(p'P~j^ip ++ {M - ui - uo 2 - E' + u[ + uo' 2 ) 



+V' + P t^- + {M + u 1 +uj , 1 - E') - r -V jjr- (M -ut-uj^- E') 

^ r - '( V -ooi- w' x - E') + r' r :^- r - (M + E') 
M M 



-V'p* -j^V (M + U1 + UJ2 — E —uj 1 — u 2 ) 



M 



0. 



(59) 



Namely, Eq. (j59j) is the representation of electromagnetic gauge invariance (the current con- 
servation, —iQ" < P'\Jfj\P >i= 0) for the instantaneous (three-dimensional) formulation. 

Now let us show a given instantaneous transition matrix element, e.g., Eq. (J48|) indeed is 
gauge invariant. 

If we multiply — iQ" to Eq.pSj) and apply : 



r d 3 q P± 




1 (2vr) 3 r< 


m y 



-iQ"Qi / T^-rr^r I — ip'p, 7/,<^p + + H>' P < 1^ P + + V ' P > TmVj 



we have 
-iQx 



fd 3 q P± j 




/ (2vr)3 Tr \ 


M L 



-^7(P- F')y + p + -y>pT{P- P')^' - p' P 7(P- />p~] •} 



(60) 
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With the properties Eqs. (|Bll IB2I IB3|) for the project operators, we may obtain the 
identities: 

r- r =A~ A = [2wiA| + J(aiM - wi + 9p)] - [2^A'+ + - ^ + g P )] 

= [-2o; 1 Ar + f («iM + wi + g P )] - [-2^A'f + § (aiS' + cu[ + q' P )\ , (61) 

A; + (/>- /")A+ = A'+^A+(M -^-uj 2 -E' + u[ + ^) , (62) 

A'+(/>- /")Ar = A; + ^Ar(M + Wl + ^ - E') , (63) 
A'f (/>- /")A+ = A'f^A+(M - a* - - Ef) , (64) 

A;-(/>- /")AI = A'f-^Ar(M + ^ + - £' - ^ - u/ 2 ) . (65) 

To apply the above identities for the project operators and with the definitions Eqs. (J23 HBl 
EUESI), one can turn the equation Eg. ((00} to the equation Eq. lJHnj) . the terms proportional 
to Q\ accordingly. For the terms proportional to Q2 it is the same. This means the gauge 
invariance (current conservation) is satisfied, and also show that it does not matter to 
multiply —iQ^ to Eq. (|48j) first and then to make the contour integration, or to reverse the 
order of the two steps, even if this two ways to do the contour integration may meet different 
poles. 



IV. DISCUSSIONS AND OUTLOOK 



As shown in the derivations, it is not so straightforward to guarantee the gauge invariance 
for the instantaneous matrix elements, so one should be full of confidence to apply them 
to various processes. We derive the Eqs. pHl |4~§j) without approximation and to start with 
a full gauge conserved four-dimensional formulation on BS wave functions, no matter the 
BS kernel being instantaneous or not (just carry out the contour integrations when the BS 
kernel is instantaneous), so it is understandable that the final instantaneous matrix elements 
Eqs. (148 1 l4*H|) are gauge invariant (the matrix element for current is conserved). Since a 
physical process, no matter it contains bound state(s) or not, must be gauge invariant, so 
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from the equations here one may conclude that the components correspond to the negative 
energy projection in the BS wave function should be also treated carefully, even the BS 
equation with an instantaneous interaction, no matter one wishes full relativistic effects are 
taken into account or not. If only the positive energy component, although being the biggest 
one in wave function, is kept, one can not have a controlled result, because it violates gauge 
invariance and most of the processes are sensitive to the gauge choice. 

If in addition to the 'positive energy' component equation Eq. (j29|) which Salpeter kept 
only, the 'negative energy' component equation Eq.()30p and the 'constraint equations' 
Eq. ljHlj) between the positive and negative energy ones for the BS wave function, as em- 
phasized here, all are applied simultaneously, then the final equations corresponding to the 
scalar components of the full BS wave function become a group of consistent equations being 
coupled each other. The fact is shown by the example Eqs. (JHH1 clearly. 

If keeping all the components of the BS equation, including the 'negative energy' com- 
ponents Eq. (j30|) with the 'constraint' Eq. (j31)) . then it is indicated that accordingly the 
effective theories, such as NRQED |7] and NRQCD [8] etc, should be affected in certain de- 
gree. It is because that the effective theories, such as NRQED and NRQCD, are established 
based on the considerations about the symmetries, gauge invariance etc, but the components 
(sectors) corresponding to the fermion and anti-fermion field are decoupled completely at 
beginning^, 0. Indeed the effective theories can be expended in 0(v) (where v = L^J i s the 
relative velocity) accurately and the gauge invariance may be kept precisely for each sector 
itself of the fermion and the anti-fermion respectively, but they can never include any of the 
special 'relativistic effects' which correspond to the relations and/or constraints between the 
fermion and anti-fermion sectors i.e. only the 'positive energy' component equation Eq. (J29|) 
itself can involve some of them, but cannot involve those which come from the 'negative 
energy' component equation Eq.fpJj) and the constraint equation Eq. tpTTI) . 

For relativistic quantum field theories, the fermion and its anti-particle are related defi- 
nitely. Dirac spinor (four components) relates the fermion and its anti-particle sector prop- 
erly but Pauli one cannot. NRQCD and NRQED are established at beginning on the Pauli 
spinors for the heavy quark sectors with fermion and anti-fermion sectors completely be- 
ing decoupled, so they cannot involve all of the connections from the relations and and 
constraints such as Eqs (pM EITj) besides Eq.(j2HJ). Note that since here it is the 'specially' 
case with instantaneous interaction shown, so the precise relations and/or constraints are 
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determined by Eqs.(j2H ISHJ) ? EU whereas in general cases they should be restricted not as 
here with the 'instantaneous interaction. If one would not ignore the negative energy com- 
ponents and wish to involve the effects in non-relativistic effective theories on corresponding 
underlying ones such as the BS equations i.e. on relativistic quantum field theories for the 
bound states, the effective theories should have certain corrections to dictate the relations 
and/or constraints between the sectors of the fermion and its anti-particle. Namely even 
each of the sectors for the fermion and its anti-particle itself is a non-relativistic with quite 
accurate corrections in 0(v) and gauge invariant, but certain relations and or constrains 
between the fermion and its anti-particle sectors should be involved in correctly, especially, 
when the relations and/or constraints are put on, then the gauge invariance should be well- 
guaranteed, as long as one does consider the fermion and anti-fermion to form bound states 
and wish to reach to a high accuracy for transition calculations between the bound states. 
As indicated by Eqs.(j2U E3 EH ESj) (although they are obtained when the binding inter- 
action is instantaneousO, this kind corrections focused here for the effective theories vary 
with the quantum numbers, such as the quantum numbers J PC (when the bound states are 
formed by a fermion and its anti-particle exactly) of the bound states. We will discuss this 
topic more precisely elsewhere 

We also think that the instantaneous formulations obtained in the paper for the radia- 
tive transitions and the BS equation should have various tests, e.g., to solve the coupling 
equations obtained from the BS equations with an instantaneous kernel, and to apply the 
obtained solutions (the spectrum and the corresponding wave functions) to the radiative 
transitions etc. Such applications of the formulations to positronium and charmonium as 
precise examples are in progress, and preliminary numerical results have been obtained 
already jl^. 
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APPENDIX A: THE GENERAL FORMULATION FOR THE 'INSTANTA- 
NEOUS' BS EQUATION 

Generally, a non-relativistic binding system can form several (even infinity) bound states, 
i.e., the BS equation have several (even infinity) solutions to correspond the spectrum of the 
binding system. When possible transitions among them are concerned, and the system is 
'instantaneous', we need to formulate the 'instantaneous' BS equation and the corresponding 
BS wave functions not only on the solution itself but also on a different solution from it. 
Namely the initial and final bound states, all of them should be formulated in one reference 
frame, i.e., both of them are either in the reference frame of the initial state P (with BS 
wave function Xp(q)) or m that, of the final state P' (with BS wave function Xp' (?'))■ 

Considering the applications, here let us describe the precise formulation of the 'instan- 
taneous' BS equation and the relevant BS wave functions Xp^) an d Xp'W' 1 ) n °t only on 
the state P but also on another one P'. In the text of the paper, we have established the 
'instantaneous' formulation of the BS equation 'set on' the state P, which may be moving or 
rest to the reference frame. Now we need to establish the formulation of the BS equation for 
the state P', but again in the reference frame of the state P (another BS equation solution 
for the bound state P'). Thus we need to divide the momentum P' and q', as done for the 
relative momentum q for Xp(g M ), into two components: a parallel component (P' p », q' p ») and 
an orthogonal one (P' p ±, q' p ±) to the momentum P (the momentum of another bound state) 
respectively. For convenience, we write the divisions for P' and q' as well below: 

P\\ ' P-L ' 

<? — Q P \\ + Q P ±, 

where P'» A = (P-P'/M 2 )P^, P'» pi _ = P'^-P'^ ph q'^ = (P-q'/M 2 )P^ and q'^ ± = q'^-q^y 
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Correspondingly, we have the Lorentz invariant variables: 

p. p> 



P' 



q' P 



M ' 

P-q' 

M P 



P' 



q'p T 



p'2 p _ p/2 



q' 2 P — q' 2 



-P 



12 



-A 



The 'instantaneous' interaction kernel is the same as Eq. (jl8|) : 

V{P^ P, q'n = V(\e p± - ggj) = V(P"*, F, q'n 
Thus the relevant 'instantaneous' BS wave function and equation are 

and integration 



r}'(q 



p-L 



kp T dkp T ds 
(2vr) 2 



v{\K L - q '^\)vAK 



(Al) 



(A2) 



(A3) 



then the BS equation becomes 

Xp>(C p q' P \) = Sf\p'SW(e±)Sf(-p' 2 n , (A4) 

where S ( j L \p' 1 fl ) and S^\~ p' 2 ^) are the propagators of the quark and anti-quark respectively 
with the momenta Eq.®. The propagators Sf(pi), Sf\— p' 2 ) and relevant projection op- 
erators of the bound state P' but in the frame of P can be written as follows. 

A' + ,p(e) A'-ip(^) 



iJSf{Jp"*^ 



+ 



Jq'p + O-iP'p — u)' iP + ie Jq' P + a^P'p + u' iP — it 



(A5) 



with 



ip 



m? - (q' p± + «;P P± ) 2 , J 2P = fitf - (q' p± - a' 2 P P± ) 2 



1 


V 




_M 


1 


'V 


2(j' 2P 


M 



-J 2p T ( m ' 2 + / -a' 2 P' 



Here A^ P (g P± ) satisfies the relations: 



(A6) 



AS(5J!l)^(^ ± ) = 



Now tppt{q'p±) are defined as: 



M 



(A7) 



(A8) 
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then 

Integrating out q' P (contour integration) on both sides of Eq. (jA4j) we obtain 

(P' P -u' lP -u' 2P ) (P P + u' 1P + u J ' 2P ) ' (A9) 

and applying the P-projection operators A' ± ip{q'^ 1 ) to Eq. (|A9J) further properly, we obtain 
the coupled equations: 

(P' P - J XP - 4p)#(0 = A'V(^V(#JA'M^) , (A10) 

(Pp + J 1P + u> 2P )ip P 7{q^) = -^ip(€±W(Q P \)A'~2p(q^) , (All) 

^-(9?J = ^ + (8?±) = 0, (A12) 
the similarities of the equations Eqs.jjZH EH EI}- 



APPENDIX B: USEFUL FORMULAE FOR THE PROJECT OPERATORS 

In this appendix, we present the formulae for the project operators which are used in the 
proving the gauge invariance for the instantaneous formulation of the radiative transitions. 

When the photon is emitted by the 1-quark Q = pi—p[, i.e. the 2-quark is kept unchanged 
(Fig.l) 5 3 (— q P ± — a 2 P'p± + q'p±), then we have u' 2 = uj 2 and A^ = Af. Furthermore it is 
easy to prove the identities for the project operators A'^p(q P± ): 

Af (9?JAf (g?J = ±^K ± (q p \) , A±(<£JA±(<£J = ±^A±(^) ; (Bl) 
A^tfJArfcft.) = (/p± + ? rP±) K( q ' P \) , Af (g^)Ar(g^) = ^Ar(g^) ; (B2) 



ATCtfJAf = Ar(g^ ± ) (/p± + ? ^ p±) , A^Aft^) = A^)^ . (B3) 

UJi U)\ 

When the photon is emitted by the 2-quark i.e. the 1-quark is kept unchanged (Fig.2), 
similar relations may be obtained. Namely the photon is emitted by the anti-quark, i.e. 
there is the condition 5 3 (q P ± — a\P'p± — q' P ±), then we have u[ = ui, A'f 1 = Af and the 
A' 2P (qp ± ), then the identities: 

A^JA^gJJ = T^A^(^) , Af(q» ± )Ai(q» ± ) = T—Af(q" P± ) ; (B4) 

Ld 2 UJ 2 
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(B5) 
(B6) 
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